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We present a model space particle-hole Green's function calculation for the quadrupole excitations 
of cold Fermi gas near Feshbach resonance using a simple model where atoms are confined in 
a harmonic oscillator potential. Both the Tamm-Dancoff and random phase approximations are 
employed. By summing up exactly the ladder diagrams between a pair of interacting atoms to 
all orders, we obtain a renormalized atomic interaction which has well defined and identical limits 
as the scattering length tends to ±cx). The experimentally observed abrupt rise in the excitation 
spectrum and its associated large decay width are satisfactorily reproduced by our calculation. 

PACS numbers: pacs 



Introduction: The Feshbach resonance was originally 
used in nuclear physics 01 , to describe a resonance in neu- 
tron scattering due to the formation of a quasi-bound 
state between the scatterer and the incident particle. 
During the past several years, this resonance was ex- 
perimentally realized in ultracold alkali gases 0. 0] ■ It 
has many interesting features and has opened a new era 
to the study of ultracold gases. It serves as an experi- 
mental control mechanism to tune the scattering length 
over a wide range, essentially from positive to negative 
infinity. Through such process, the inter-atomic effective 
interaction could also be experimentally controlled. As 
a result, ultracold alkali gases have provided us a unique 
many-body system with tunable interactions, and they 
have been the focus of many theoretical studies (See, e.g. 
SB)- As the magnetic field is tuned across the Feshbach 
resonance, ultracold Fermi gases have been observed to 
undergo the well known BCS-BEC crossover 0, 0, II] ■ 

Collective oscillations of trapped cold Fermi gas have 
been recently measured |^ Il0j| , with important and very 
interesting findings for the coUisional-, coUisionless- and 
particularly the transitional-regime in between. In the 
coUisional regime, the oscillation frequency is lower than 
the trap frequency and it varies continuously across the 
Feshbach resonance where the scattering length tends 
to ±oo. In the coUisionless regime, the frequency is es- 
sentially constant, close to the trap frequency. In the 
transitional regime, however, an abrupt rise of frequency 
together with a strong damping rate were observed. To 
explain these findings has been a challenging task. In 
the coUisional regime, pioneering hydrodynamical stud- 
ies 0,^3 have been carried out. 

In this Letter, we would like to carry out a micro- 
scopic study for the quadrupole excitations of trapped 
cold Fermi gas, aiming primarily at the transitional and 
coUisionless regimes. We shall use a particle-hole Green's 
function framework which has been commonly used in 
nuclear physics 0, 113. ITsI l . When applying this method 
to nuclear systems, nucleons are placed in an oscillator 
well and particle-hole excitations are treated to all orders 
using either the Tamm-Dancoff or random-phase approx- 
imations, commonly referred to as TDA and RPA respec- 



tively. We shall use a similar approach for the above gas 
system, with its atoms trapped in a harmonic oscillator 
potential. 

In the following section, we shall first briefly describe 
the model space particle-hole Green's function method 
[l^ which we shall employ. A specific feature of the 
atomic interaction for cold Fermi gases is that it is 
divergent at Feshbach resonance, [j, E, 0, 0, 
This interaction is not suitable for being used directly 
in microscopic calculations; it needs to be renormal- 
ized (or tamed) beforehand. We shall discuss that the 
well known Brueckner reaction matrix method, which 
has been widely and successfully used in nuclear matter 
[23L 0, and finite nuclei , is very useful in this re- 
gard. We shall use this method to derive a renormalized 
atomic interaction which is smooth and well behaved at 
the Feshbach resonance. Our results together with some 
calculational details will be reported and discussed in the 
last section. A summary will also be presented. 

Model space particle-hole Green's function 
method: In this section we briefly describe the model 
space particle-hole {ph) Green's function method 0| for 
calculating the quadrupole excitations of cold Fermi gas. 
We consider atoms as spin 1/2 fermions confined in a 
spherical harmonic oscillator potential of oscillator spac- 
ing hujosc and oscillator length aosc = y^h/muJosc- In our 
present calculation, we assume the atoms remain normal 
throughout. As pointed out in trapped fermi gas 

experiments, the superfluid component always resides in 
the interior of the trap and had little contribution to the 
quadrupole excitation. We shall use hujosc and aosc re- 
spectively as the units for energy and length. 

We treat the trapped g closed shell system of Np 

closed shells. For example, the Np — 6 system has 112 
atoms and in its unperturbed ground state the lowest 
6 oscillator shells are completely filled. The excitation 
energies, denoted by i?„, of th e sy stem are calculated by 
solving the ph RPA equation [ij) UM 

AXn+BYn = EnXn, 

-B*Xn-A*Y,, = EnYn. (1) 
Here A represents the sum of the unperturbed ph energy 
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gap and the TDA vertex function, namely 

A{ph,p'h') = iep + Sip)-eh-Sih))6ph,p'h'+^iph,p'h'), 

(2) 

where S denotes the one-body vertex function and E the 
two-body ph vertex function. They are composed of ir- 
reducible diagrams as illustrated in Fig. 1. Here each 
dashed-line represents the vertex of the effective atomic 
interaction V. Note that the diagrams are time ordered 
with retarded single particle propagators Q . B denotes 
the ground-state correlation vertex, whose diagrammatic 
structure is also illustrated in Fig. 1. All the diagrams 
contained in the vertex function must be irreducible in 
the sense that any intermediate state of such diagrams 
must be orthogonal to the ph model space . 






FIG. 1: Diagrams for one- and two-body vertex functions. 
The dotted line vertex represents a V interaction. 

The transition amplitudes X and Y of Eq.(l) are 
X„{ph) = (*„|ata,,|^'o) and Y„{ph) = (^-^Jat ap|*o), 
where "^1^ and ^'g represent the nth excited- and ground- 
state wave functions respectively. When the ground state 
correlation vertex B is suppressed, Eq.(l) becomes the 
TDA equation. We consider that the atomic gas is 
trapped by a harmonic oscillator potential. Thus e of 
Eq.(2) is the harmonic oscillator single particle energy. 

In carrying out the RPA/TDA calculations, a re- 
stricted model space near the Fermi surface is employed. 
We use a hole space consisted of three major shells im- 
mediately below the Fermi surface, namely shells Np to 
{Np — 2), and for particles one major shell immediately 
above. The external lines p-,h,p' and h' of the vertex 
function diagrams of Fig. 1 are all confined in this space. 
That we use a restricted model space requires a model- 
space renormalization, so that the effects from the or- 
bitals outside the model space are taken care of. The 
inclusion of the higher order diagrams in the vertex func- 
tions is for this renormalization purpose. As an example, 
let us consider the third order diagrams of Fig. 1. They 
all have a common "ladder" structure, with repeated in- 
teractions between a pair particle lines. The intermediate 
states of them include all the particle states outside the 
model space, so that our ph calculation within a small 
model space have renormalization contributions from the 
space outside. As discussed later, this type of "ladder" 



diagrams can be summed up to all orders, giving rise to 
the reaction matrix interaction. 

Separable reaction matrix interaction: In car- 
rying out our microscopic particle- hole calculations, we 
must first have an atomic interaction on which our many- 
body calculation is based. This interaction [l. ITgL ItR lia . 
Egf has a specific form and is commonly given as 

V = 47r—aJ{r), (3) 

TO 

where m denotes the atomic mass, and as is the s-wave 
scattering length which can be tuned across the Fesh- 
bach resonance, essentially from ~oo to -foo. The above 
is an effective interaction for one hyperfine channel, in 
the sense that we consider atoms as of only one hyper- 
fine state with the other hyperfine states having been 
integrated out by renormalization. In our calculation we 
assume atoms interacting with s-wave interactions only. 

The above interaction is clearly divergent at the Fesh- 
bach resonance (the unitary limit) where —^ ±oo, and 
there is difficulty in using it in microscopic calculations. 
As described in section 2, to carry out the TDA/RPA cal- 
culations, we need first calculate the vertex functions S, 
S and B. As illustrated in Fig. 1, these vertex functions 
are composed of irreducible diagrams. Near the Feshbach 
resonance, the above atomic interaction is divergent and 
consequently every diagram in the figure is divergent; cal- 
culation can not go on unless we use a different approach. 
This type of difficulty has been well known in nuclear 
matter theory and has been overcome by the Brueckner 
reaction matrix theory [23L 123 . l25l | . There each vertex of 
the bare nucleon-nucleon interaction Vnn is nearly diver- 
gent because of its very strong repulsive core. However, 
the Brueckner reaction matrix G given by the all order 
summation of the Vnn interactions between a pair of 
nucleons is finite and well behaved. 

We apply a similar reaction matrix approach to our 
present atomic gas calculation. Let us consider the di- 
agrams for the one-body vertex function S of Fig. I. 
Clearly they have the structure of a geometric series and 
can be summed up to all orders. Let us define a reaction 
matrix R as 

Q Q Q Q , . 

R = V + V-V + V-V-V + ... = V + V-R, (4) 

e e e e 

where e stands for {uj — Hq + zO+) and Q the Pauli exclu- 
sion operator. (Recall that we use time-ordered diagrams 
with retarded single particle propagators.) Hq is the un- 
perturbed Hamiltonian for the pair of interacting atoms, 
and uj is the starting energy which we shall discuss later. 
The role of Q is to ensure that the intermediate states of 
R must be outside of the chosen model space. In terms of 
R, the sum of all the ladder- type diagrams as shown by 
S{p) of Fig.l is now {ph\R\ph). Similarly the sum for T, 
is {ph'\R\p'h) and for B is {pp'\R\hh'). We shall include 
only diagrams first order in R for the vertex function as 
listed in Fig. 2. 
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FIG. 2: First-order R-matrix (wavy-line vertex) diagrams for 
one- and two-body vertex functions. 



In calculating R, we shall use an angle average approx- 
imation [24J for the Pauli operator, namely Q is approx- 
imated by Qaviq, kp) where q is the relative momentum 
and kp is an average Fermi momentum which can be es- 
timated by using for example a local density approxima- 
tion. In addition, we shall use plane-wave intermediate 
states and an average starting energy lo. Accurate meth- 
ods for calculating the nuclear reaction matrix have been 
developed , and we plan to use them in a future work. 
With these approximations, the above R matrix becomes 



R(k,k',uj) = V{k,k' 



V{k,q) 



uj - q^ + iO+ 



R{q,k',uj)q^dq, (5) 



where k,k' ,q are the relative momenta. 

Let us now consider the above R matrix for the scat- 
tering length dependent potential of Eq.(3). Schemati- 
cally, R is given hy R = V/{1 — VQ/e) and may have 
well defined limits as V ±00. But how to obtain 
them accurately is a difficult task in general, and it may 
not be feasible to obtain them reliably using numerical 
methods. We feel that R has to be solved analytically in 
order to obtain the above limits reliably. In this regard, 
we resort to the separable potential approach which has 
been widely used in nuclear physics 20, 21, 22]; the reac- 
tion matrix with this approach can be analytically solved. 
We shall use a separable potential which corresponds to a 
simulation of the potential of Eq.(3) by a short range non- 
local Yukawa potential of the form ase~^^e~^^ /(/i^rr'), 
being the range of the interaction and r and r' the 
inter-atomic radial distance. We shall use large fi to sim- 
ulate short range potential. In momentum space this 
potential becomes 



Vik,k') = a,af{k),f{k')- f{q) = a - 8^-^' 



/i + q-^ 

(6) 

With this separable potential, the reaction matrix, 
which is complex, is also separable and is obtained ana- 
lytically as 



i?(fc, k', uj) = [^^re{^^, u:) + iT1„^{^l, u;)]f{k)f{k'), (7) 



r]re{p-,U!) = 



Difi^io) = 1 + 
X{fi,uj) 



(8) 



1 — agal{fj,, Lu) 

Note that rjim = if w < 0. In the above, / is the integral 
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Qav{k,kp) 



{uj~k^){fi'^ + k'^y 



rk'dk. 



(9) 



where P denotes principal value integration. In a recent 
study of fermionic atoms in optical lattices |23| , a Gaus- 
sian separable interaction is employed and its strength 
has similar dependence on as the above A. 

It is readily seen that R has well defined limits as — s- 
±00. Only the factor A of Eq.(8) is dependent on a^, and 
at the above limits it becomes 



A±oo(/^,t^) = 



-1 



I{^l,uJ)' 



(10) 



Thus our i?-matrix has a common well defined limit at 
the Feshbach resonance, being the same whether as — > 
-l-oo or as —00. 

For dilute systems, kp should be small. We have cal- 
culated the above R matrix for several values of kp, and 
found that the results for R calculated with kp ranging 
from to 1 are quite close to each other. For the case of 
kp = Q and uj > Q, the integral of Eq.(9) becomes 



I{n,uj) 



UJ — jl 



(11) 



which is only weakly energy dependent if lo is much 
smaller than /i^ . Thus our _R- matrix interaction is nearly 
energy independent for cold (small uj) dilute (small kp) 
system with a short-range (large fi) interaction. 

Results and discussion: We now describe some de- 
tails about our calculations for the energy and decay 
width of the quadrupole excitation. They are obtained 
from the solutions of the RPA/TDA equation (1), us- 
ing the vertex functions given by the diagrams of Fig. 
2. These diagrams are calculated from the _R-matrix in- 
teraction, using the methods detailed in 26]. Since R 
is complex, Eq.(l) is now a complex equation and its 
eigenvalues £"„ are generally complex. We write En as 
{ReEn — iTn), with the decay width given by r„. In Fig. 3 
we present results for the lowest 2+ state of an Np—8 
(240 atoms) model system, using parameters kp—0.3, 
uj—1 and /i=20. (We have also calculated Np=6 and 7 
systems, with results both being highly similar to those 
of Np = 8.) These parameters should be suitable for 
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FIG. 3: Energy and decay width of the quadrupole excitation. 

cold dilute systems interacting with a short range inter- 
action. The above choice has the advantage of making 
i?-matrix weakly energy dependent as we have fj? >> to. 
As shown, our calculated energy at Feshbach resonance 
(FR) is about 1.7 huosc with nearly zero width, both 
being smooth functions there. 

The shapes of our calculated spectrum and decay 
width are worth noting. Starting from the left end of Fig. 
3, ReE exhibits a gradual drop followed by an abrupt rise 
and finally remains nearly constant. At both ends, ReE 
approaches to 2. Most importantly, the decay width F 
is peaked concurrently with the abrupt rise, both hap- 
pening right at a\ . These specific features of our results 



agree remarkably well with experiments 0, . 

The above abrupt rising feature is mainly due to the 
i?-matrix interaction. The factor A of Eq. (8) has a pole 
when Os equals to ax — —l/{al). It is readily seen that 
at a\ the imaginary part of the R matrix is peaked with a 
Lorentzian form, while its real part vanishes. In addition, 
R is positive to the left of {—ciosc/ o,\)^ and negative on the 
other side. These properties of R determine the general 
shapes of our results. For small kp and lo, a\ is close to 
aosc/47r. In fact for kp=0 and uj ~ 0, ax is exactly equal 
to aosc/47r, independent of /j. As shown in the figure, the 
abrupt rising takes place very close to 47r. We note that 
RPA and TDA give nearly identical results far from ax- 
Near ax, ReE for RPA tends to negative indicating the 
system being unstable. 

Last but not least, F is discontinuous at ax as shown. 
We have found that the wave functions X and Y of Eq. 
(1) are also discontinuous at ax- We plan to further study 
the nature of this "transition" in a future work. 

In summary, we have carried out a particle-hole 
shell model calculation for the quadrupole excitations 
of trapped cold Fermi gas. An essential step is the use 
of a reaction matrix renormalized interaction. Using a 
separable potential approach, this interaction has been 
derived analytically and it is a smooth function at Fes- 
hbach resonance although the atomic interaction before 
renormalization is divergent there. Our results demon- 
strate that, in between the coUisional and coUisionless 
regimes, trapped cold Fermi gas undergoes a transition 
with abrupt variations in both energy (frequency) and 
decay width (damping rate), in good qualitative agree- 
ment with experiments. 
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